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ABSTRACT 

In this paper we present a method to analyse absorption line spectra of a galaxy designed 
to determine the stellar dynamics and the stellar populations by a direct fit to the spectra. This 
paper is the first one to report on the application of the method to data. The modelling results 
in the knowledge of distribution functions that are sums of basis functions. The practical 
implementation of the method is discussed and a new type of basis functions is introduced. 

With this method, a dynamical model for NGC 3258 is constructed. This galaxy can be 
successfully modelled with a potential containing 30% dark matter within lr e with a mass 
of 1.6 x 1O 11 M . The total mass within 2r e is estimated as 5 x lO n M0, containing 63% 
dark matter. The model is isotropic in the centre, is radially anisotropic between 0.2 and 2 
kpc (0.88r e ) and becomes tangentially anisotropic further on. The photometry reveals the 
presence of a dust disk near the centre. 

Key words: methods: numerical - methods: statistical - galaxies: kinematics and dynamics - 
galaxies: elliptical and lenticular, cD - galaxies: individual (NGC 3258) - galaxies: structure 



1 INTRODUCTION 

It is well known that in galaxies stellar populations and stellar dy- 
namics cannot be studied completely independently; the determi- 
nation of the line profiles and kinematic parameters depends on 
assumptions about the dominant stellar population, while the de- 
termination of the chemical properties depends on the absorption 
strenght of the spectral lines. 

In general, stellar kinematics are retrieved by comparing 
the observed spectra with the convolution of a stellar tem- 
plate spectrum with a model line-of-sight velocity profile, see 
iDe Bruvne et alj 1200 3l) for a discussion of some standard tech- 
niques. The stellar template belongs to the same spectral class as 
the dominant stellar population in the galaxy, or can be a mix of 
different stellar type spectra. Often, a best matching stellar mix is 
determined for the spectrum at the centre, and this mix is then used 
throughout the whole galaxy. In som e cases, the stellar mix is al- 
lowed to vary with radius, (e.g. lCarter et aljll999l) ). However, the 
only criterion used for obtaining a proper mix is a good fit, without 
any physical considerations. Therefore, it is acknowledged that this 
procedure does generally not yield a reliable population synthesis. 

The kinematics obtained from a parametrization of the line- 
of-sight velocity distributions can be used as input for a dynamical 
model. Within this multistep scheme of dynamical modelling, it is 
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acknowledged that dynamical evidence for dark matter in ellipti- 
cal galaxies requires the use of more information contained in the 
line-of-sight velocity distribution (hereafter LOSVD) than simply 
the mean velocity and velocity dispersion. The reason for this is the 
existence of a so called mass-anisotropy degeneracy. In a dynami- 
cal model it is possible to mimic the influence of a dark matter halo 
by tangential anisotropy. Therefore, also the anisotropy informa- 
tion in LOSVDs should be used. However, there is no uniformity in 
the parametrization of these LOSVDs and different algorithms may 
well have different biases lloseph et all200lh . Moreover, it is not 
always straightforward to include these parameters in a fit, since not 
all of them are linearly dependent on the distribution function (here- 
after DF) iGerhard et al]il998t) . lDe Bruvne et ai]<200ll) 1. Another 
option is to use the complete LOSVD as dynamical input. This 
omits the need to approximate the LOSVD with a number of pa- 
rameters, but in combination with the Schwarzschild orbit superpo- 
sition method, this requires a careful smoothing both of the LOSVD 
and the dynamical model iGebhardt et al. 2000). But there is still 
another way of dealing with the kinematic information of a galaxy, 
through direct modelling of the observed spectra. 

In this paper we adopt a method to analyse spectra from ellip- 
tical galaxies that explores the dynamics of the galaxy and simulta- 
neously offers a way to study the stellar populations in that galaxy. 
The method originates from the field of dynamical modelling and 
is outlined bv lDe Riicke & Deionghel <1998l) . hereafter DD98. The 
main idea is that the different stellar populations that contribute to 
the integrated galaxy spectrum do not necessarily share the same 
kinematic characteristics. Hence, in terms of dynamical modelling, 
they should be associated with different distribution functions. A 
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simultaneous determination of the best stellar mix and the best dis- 
tribution functions associated with them de facto amounts to a dy- 
namical model and a stellar population synthesis. 

This paper is the first to present the application of the mod- 
elling method to data. The aim is to determine the gravitational po- 
tential of NGC 3258, by modelling the galaxy spectra using differ- 
ent stellar templates. The modelling strategy is outlined and com- 
mented, and a new family of dynamical components is presented. 

This is the first paper in a series of two, where a deeper in- 
terplay between dynamics and population synthesis is explored. In 
this first paper, the attention will be focused on the dynamical mod- 
elling aspect, while the second paper will go deeper into the stellar 
population analysis. 

In the next section, the observations and data reduction steps 
are presented. The modelling strategy is outlined in section 3. Sec- 
tion 4 presents the dynamical components that were used. In sec- 
tion 5 the results of the modelling are presented. Section 6 contains 
a discussion of the results and the conclusions can be found in sec- 
tion 7. 



2 OBSERVATIONS AND DATA REDUCTION 

The El galaxy NGC 3258 is one of the brightest members of 
the Antlia Group. In its neighborhood, NGC 3257 can be found 
at 4.5 arcmin (about 21r e ) and NGC 3260 is at 2.6 arcmin 
(about 12r e ). We measured a systematic velocity of 2709 km/s, 
this puts NGC 3258 at a distance of 36.12 Mpc (taking H = 
75km/s/Mpc). This means that a linear size of 1 kpc at the dis- 
tance of NGC 3258 translates into an angular size of 5.7". 

NGC 3258 was observed with the ESO-NTT telescope in the 
nights of 27-28/2/2000. Spectra of the major axis covering the Ca 
II triplet around 8600 A were taken. During the observations in 
2000, also imaging of NGC 3258 in the R band was performed. 
The detector used for this observation run was a Tektronix CCD 
with 2048 x 2047 pixels, 24pmx24^m in size and with a pixel 
scale of 0.27"/pix. Another observation run, in 1999, was used to 
obtain images in the B band. 

2.1 Photometry 

R band observations (in total 600 sec integration time) of 
NGC 3258 were obtained at the ESO-NTT using the red arm of 
EMMI during the spectroscopic run in februari 2000. An addi- 
tional B band image (600 sec integration time in total) was obtained 
in february 1999 in the same configuration. The detector used for 
the B band images was a Tektronix CCD with 1024x 1024 pixels, 
24/xmx24/im in size and with a pixel scale of 0.37"/pix. 

The seeing (judged from a number of stars in the images) was 
1.1" for the R band images and 1.7" for the B band images. The 
photometric zero point was established by comparison with pub- 
lished photometric data available on NED. Figure \\\ shows an im- 
age of NGC 3258, isophotes are plotted in figure|2| In this figure, 
the central isophote seems to be a little displaced compared to the 
other isophotes. Both B and R images show a slightly off-centre 
nucleus, with a shift of about 0.45" (or 0.07 kpc) in the horizon- 
tal direction and a shift of about 0.1" (or 0.01 kpc) in the vertical 
direction. This is illustrated in the upper panel figure [5] where the 
shift of the centres of the isophotes is shown. 

The photometric parameters were derived by fitting elliptic 
isophotes to the image. The intensity along the isophotes was ex- 
panded in a fourth order Fourier series, in order to reveal deviations 
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Figure 1. B band image of NGC 3258, the size of the region is 110" x 110' 
North is up and east is to the left. 
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Figure 2. Isophotes of NGC 3258. 
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Figure 4. Inverted image of the central part of the unsharp masked R image. 
The light grey band in the centre is a central dust disk. 



Figure 3. Photometric parameters for NGC 3258. Left column, from top to 
bottom: surface brightness, position angle, ellipticity and B-R. Right col- 
umn, the coefficients indicating deviations from pure elliptic isophotes. 



off-centre nucleus appearing in the B and R images is probably due 
to the obscuring effect of the dust lane. 



from pure elliptic isophotes, like in e.g. iDe Riicke et alJ 120031) . 
In figure [3] the main photometric parameters are shown: surface 
brightness, position angle, ellipticity and B-R profile as derived 
from the B and R surface brightness profiles in the left column. 
The coefficients indicating deviations from pure elliptic isophotes 
are shown in the right column. The values for these coefficients in- 
dicate that NGC 3258 is an almost perfect elliptical galaxy, with an 
ellipticity between 0. 1 and 0.2. 

A Sersic profile (I(r) = Io exp[— (r/ro) 1 ''™], with 7(r) the 
surface brightness profile, Io the central surface brightness, rn the 
scale length, an n the Sersic parameter ICiotti et alll995h was fit- 
ted to the surface brightness, the result of the fit is overplotted (in 
solid line) in the upper left panel of figure|3| The observed surface 
brightness profile can be well matched with a Sersic profile with 
Sersic parameter n = 5.47 ± 0.07. The effective radius is 12.95" 
or 2.26 kpc. 

Values for B-R (derived from the photometric profiles) are 
shown in the lower left panel of figure [3] Within the central 5", 
there seems to be a slight reddening of the nucleus. B-R changes 
from about 1 in the most central measurement to about 0.85 at 5". 



2.2 Central dust component 

Figure [4] shows the central part of the unsharp masked R image 
with a clearly visible dust feature. The estimated radius of the ma- 
jor axis of the disklike feature is about 1" - 1.5" (0.17 kpc - 0.26 
kpc). Figure [5] shows a contour map of this region (lower panel) 
and a cut along a horizontal line through the centre (upper panel). 
This line does not coincide with any of the photometric axes. The 
effect of the dust lane on the apparent brightness is clearly visible 
in the upper panel of that figure. The brightest central point does 
not coincide with the centre of the isophotes for the main body of 
the galaxy. The same behaviour is seen in the B band images. The 



2.3 Spectroscopy 

For the spectra, grating #7 was used, having a dispersion of 0.66 
A/pix. A slit width of 1.5" yielded a spectral resolution of 3.67 A 
FWHM, resulting in an instrumental dispersion of about 54 km/s in 
the region of the Ca II triplet. For NGC 3258, several exposures of 
3600 sec were taken (in total 1 1 hours). 

A number of standard stars (G dwarfs and K and M giants) 
were also observed, they are listed in tableQ 

Standard reduction steps were applied to these spectra with 
ESO-MIDAS \ Dedicated exposures, taken during both runs, were 
used for bias correction, dark correction and flatfielding. Cosmic 
ray hits were removed using a top hat filter. After this correction, 
remaining cosmics were also removed by hand. 

For the wavelength calibration, lamp spectra were taken just 
before or after each of the spectroscopic observations of the galaxy 
or the stellar template stars. For each row in the spectra a polyno- 
mial was constructed to transform the pixel scale into wavelength 
scale. The calibrated images with the Ca II triplet have a step of 
0.3A. 

Airmass correction is applied, using the mean value of the air- 
mass at the beginning and end of the exposure. The contribution 
of the sky to the spectra was estimated from an upper and lower 
region of the image, where there was no contribution of the light 
of the galaxy or template star. Several galaxy spectra were reduced 
separately, aligned and combined into one galaxy spectral image. 

2.4 Kinematic parameters 

In figure [6] kinematic parameters are presented. These parameters 
are derived from the spectra, using a x 2 minimization technique in 
order to find the parameters of a truncated Gauss-Hermite series 

1 ESO-MIDAS is developed and maintained by the European Southern Ob- 
servatory 
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Figure 6. Kinematics obtained from the spectral feature around 8600 A. 



HD 21019 


G2V 


HD 115617 


G5V 


HD 102070 


G8III 


HD 117818 


KOIII 


HD 114038 


K1III 


HD 44951 


K3III 


HD 29065 


K4III 


HD 99167 


K5III 


HD 129902 


M1III 


HD 93655 


M2III 



Table 1. Observed stellar templates. 



Figure 5. Upper panel: displacement in x (A(x c )) and y (A(y c )) direction 
of the centre of the isophotes as function of the major axis of the isophotes. 
Middle panel: The intensity profile along the cut with a horizontal axis 
through the galaxy centre (dotted line in lower panel). Lower panel: con- 
tours of the central part of the unsharp masked R image. The photometric 
major and minor axis are plotted in dashed line. 



for the second order moment, (v p ) and a v are also called the 'cor- 
rected' values. 

The se data are consistent w ith the kinematic paramete rs pre - 
sented by Pelleg rini etalli 19971) and Koprolin & Zeilinger (2000). 



that gives the best approximation to the LOSVD. The spectra were 
spatially rebinned to satisfy a minimum S/N. Only data points from 
the side of the galaxy with mean streaming away from the observer 
are shown. A best fitting template mix (i.e. the one that gives the 
smallest x m the fitting) was determined for the central data point 
and this mix was then used throughout the whole galaxy. The spec- 
tra are analysed with a mix of a G5V dwarf and a K4III giant. 



method of 


van der Marel & Franx 


De Riicke et al. 


12003), first fitting 



mine (v) and a, and then keeping these lowest order parameters 
fixed to fit the other coefficients. 

If (v), a, h$ and h± are the lowest order parameters for the 
LOSVD, the true moments are approximated more accurately by 



(v p ) = (v) + V3ah 3 

for the first order moment and 

a v — <t(1 + Vo/u) 



(1) 



(2) 



2.5 Other observations 

NGC 3258 is a weak but extended radio source with a total flux 
of 18.3 mjy (VLA observations at 4.9 GHZ (6 cm) within 45") 
( Sadle r et alll989h . where the core « 5") flux is only 7% of the 
total flux. A hot gaseous corona arou nd NGC 3258 emits about 
3 x 10 41 ergss _1 ( Fab iano et all 1 992). There are no total mass es- 
tima tes based on X- r ay obs ervations available. 

iBregman et alj Q998) studied this galaxy in the far-infrared, 
and derived a dust temperature of 24.5 K and a dust mass of 
6 x 10 6 Ai© from IRAS observations at 60 and 100/im. This is not 
unlike what is found in spiral galaxies. The mid-infrared properties 
of NGC 3258, based upon ISOCAM observations at 6.75, 9.63 and 
15Mm, indicate the presence of a hot T ~ 260K) dust component 
(Ferrari et al. 2002). A mass of A&Mq of hot dust within a radius 
of 5.5 kpc was derived. This is at the lower end of the typical range 
of 10-400 Mq reported for early-type galaxies. 

Fur thermore, the galaxy shows [QUI emi ssion 
jBettoni & Busonll987l) and weak HI emission (Jen kinJl983l) . 
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3 DYNAMICAL MODELLING STRATEGY 

There is a large variety of ways to make a dynamical model and 
there also different strategi es to use thes e models. For a recent re- 
view, see Deionahe & De Bruvne 120031) . 

The data used for the modelling are the projected surface 
brightness and a set of spectra of the galaxy, taken at several posi- 
tions along the major axis. The projected surface brightnes is used 
to obtain a graviational potential, this is described in more detail in 

rrn 

An important issue for the modelling strategy applied in this 
paper, is the form of the data that are used to determine the DF. 
In this case, the fit is done directly to the spectra (after necessary 
data reduction steps). This means that the observed quantities that 
have to be reproduced by the model are pairs (wavelength, flux). 
The use of spectra as input for dynamical information is in contrast 
with more frequently used modelling techniques where a number of 
kinematic parameters are first derived from the spectra. The mod- 
els are then fitted to these parameters that are moments of the dis- 
tribution function. The modelling procedure used in this paper is 
described in detail in DD98. The input spectra are reproduced by 
a dynamical model that uses spectra of different template stars. An 
important characteristic of these models is that they provide a DF 
for each population that is used in the model. 

This is a function which gives the probability to find a star on 
an orbit characterized by the integrals of motion E (the energy), L 
(the total angular momentum) and L z (the component of the angu- 
lar momentum along the rotation axis of the system). This function 
holds all information on the dynamical state of the galaxy. 

The resulting DF for each type of star is a sum of basic compo- 
nents and is con structed with a quadratic programming algorithm 
lDeionghel l989). There are two main requirements for the DF: (1) 
the difference between the observed quantities and the quantities 
calculated from the DF should be minimized (least squares mini- 
mization) and (2) the DF should be a positive function. 

Modelling the spectra instead of kinematic parameters causes 
a considerable increase in the number of data points and hence an 
increase of the required computing time, but there are a number of 
advantages over using kinematic parameters. 

• M ore traditional modelling involves a two step process (e.g. 
iDe Bruvne et alj fcOOll) ): first determining kinematic parameters 
out of observed spectra and then calculating a model. Here this 
is turned into a one step process. Omitting the step of determin- 
ing kinematic parameters and using the galactic spectra as in- 
put has the consequence that the kinematic information used as 
input for the modelling is independent of any parameterization. 
Also in the more traditional two step modelling process, there are 
some non-parametric methods to derive kinemtical profiles, see 
IDe Bruvne et"alll2003l) . 

Of course, the dynamical profiles that are a result from the mod- 
elling depend on the assumptions concerning the gravitational po- 
tential and the dark matter halo, but this is the case for every dynam- 
ical modelling technique. In this paper, the gravitational potential is 
assumed to be spherical. This assumption is not a conceptual issue 
but is justified by the shape of the NGC 3258. 

• Furthermore, since the basis distribution functions that are 
used in the fit are chosen freely by the algorithm from a large li- 
brary, also the modelling procedure itself is non-parametric. The 
smoothness of the model is guaranteed by the use of smooth ba- 
sis functions. Hence, considering also the previous remark, this 
method can be regarded as a non-parametric modelling method. 
Moreover, the model has considerably more freedom in determin- 



ing LOSVDs than can be achieved with a limited number of param- 
eters as is the case when kinematic parameters are determined in a 
separate process. 

• The modelling process based on quadratic programming, uses 
a x 2 ( see equation (10) in DD98) that has the statistical meaning 
of a goodness of fit indicator if the noise on the spectra is Poisson 
noise and the data points are weighted accordingly. 

• The Hessian matrix (see equation (13) in DD98) involved in 
the quadratic programming can be used to calculate error bars on 
the distribution functions that have a statistical meaning and give 
an idea of the uncertainty on the model and on its moments. Tests 
with synthetic spectra (DD98) show that the error bars are a reliable 
measure for the spread between input model and fit, as caused by 
the adopted Poisson noise. Moreover, if the component library is 
sufficiently complete the systematic errors are negligible compared 
to the random errors. 

• If the dynamical modelling uses different template stars, it can 
result in a population synthesis. 

• The composition of the stellar mix is allowed to vary smoothly 
as a function of radius and is therefore directly connected to the 
internal dynamics. This minimizes stellar template mismatch, and 
allows in meantime to obtain a physical stellar population decom- 
position. 

• Information at different wavelengths can be combined in an 
elegant way. This is an advantage of using kinematic information 
on the level of spectra. It is not clear how to interpret profiles of 
the same kinematic parameters taken from different parts of the 
spectrum that give different results J van der Mare! 1994) . 

3.1 Templates 

Five template stars that have discernible spectral properties were 
used: a G2V dwarf, a G5V dwarf, a K1III giant, a K4III giant and 
a M1III giant (see tableQfor their identification). 

The dynamical models were created in an iterative process. 
First, models were calculated with a single template star, using a li- 
brary of about 100 dynamical components. Afterward, the compo- 
nents that are selected for the single template models are combined 
into component libraries that contain a sample of dynamical com- 
ponents for different stellar templates. With these mixed libraries, 
new models were calculated. 



4 DYNAMICAL COMPONENTS 

Using a spherical potential in combination with rotating compo- 
nents for the DF, it is possible to obtain a flattened luminous density 
distribution, as is required to be able to fit the two dimensional sur- 
face brightness of this El galaxy well. We consider three integral 
models: F(E, L, L z ) for which oy, cr$ and ag are different. These 
models are linear combinations of three types of components. 

4.1 Continuum terms 

The continuum of the weighted sum of template spectra does not 
necessarily fit the continuum of the galaxy spectra. Therefore, low 
order polynomials are included in the component library to com- 
pensate for the small differences between both continua. Such a 
continuum component has the form 

( In A -In Ax \° ( j R {r p ) - j R (r c ) V 
g(lDA ^H mA 2 -lnAj { JM J ' (3) 
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with a and (3 real numbers, r c the cutoff radius (the component 
is set to zero for larger radii), in(0) the central value of the sur- 
face brightness, jR(r c ) the surface brightness at the cutoff radius, 
ifl(^p) the surface brightness at the position of the line of sight, 
A the wavelength at which the spectrum is evaluated, Ai and A2 
bound the spectral region that is used in the fit. The value of a 
should be small so that these components only fit the continuum 
and not the absorption lines. 



4.2 Fricke components 

These components have the following augmented density 



p(tp,r) = (ip - 



2/3 



ip ^ Eq, elsewhere 



(4) 



with a and (3 positive integers, r c a scale length and Eq the energy 
level at which the component is truncated. All of these components 
are either isotropic (/3 = 0) or tangentially anisotropic (compo- 
nent parameter (3 > 0). However, the coefficients for these compo- 
nents can be negative, hence also radial anisotropic models can be 
constructed. These are spherical two integral (E, L) components. 
They have equal velocity dispersions in the tangential directions, 
o> = ere, but a different radial dispersion in the case of non- 
isotropic components. A more detailed discussion of this family 
of components can be found in DD98. 



4.3 Rotating components 

These components add rotation to the model. Due to the ordered 
motion of the stars around the component's symmetry-axis, the tan- 
gential dispersion is usually lower than the radial dispersion. These 
rotating components are axisymmetric two integral (E, L z i) com- 
ponents (with z' the symmetry axis of the component), for which 
or = cr z , and the dispersion in the rotation direction cr^ is differ- 
ent. 

In this section, a family of components is introduced that 
should enable one to reproduce both the rotation of the bulk of the 
galaxy and the kinematics of a small subsystem, rotating around 
an axis other than the main galaxy's rotation axis. The distribution 
function of the components of this family is given by 



F (E,L Z ,) = — (E - E ) a (i Iot L z 

for E > E 
and i rot = 1, L z , 
ori rot = -1, L z 



L f 



> and L z i > L 
< Oand \L Z ,\ > L , (5) 

with the energy per unit mass E = tp(r) — v 2 M /2 — L 2 Z , /(2-uj 2 ) 
and the z'-component of the angular momentum L z i = zuv^. The 
radial distance to the centre is measured by r, the distance to the 
centre measured in the equatorial plane is denoted by vj. The ve- 
locity component in the meridional plane is vm = \J v% + v 2 , , 
the velocity component perpendicular to the meridional plane is 
given by vj,. Here, (x' , y' , z') is a Cartesian reference frame, at- 
tached to the component, with the z'-axis the rotation axis of the 
component. The overall gravitational potential ip(r) is taken to be 
spherical symmetric although the components themselves have an 
axisymmetric mass distribution. This is acceptable since a small 
kinematically distinct core probably does not significantly influ- 
ence the global potential. Also, a slowly rotating elliptical galaxy 
can still be approximately round, justifying the use of a spherical 
potential. An axisymmetric potential generally allows E and L z 



as integrals of motion, with L z the component of the angular mo- 
mentum parallel to the potentials symmetry axis. Therefore, these 
components can also be used in an axisymmetric potential if their 
rotation axis z' coincides with the symmetry axis z of the potential. 

All the stars have the same sense of rotation, i.e. if i rot = 1 the 
stars rotate counter-clockwise, if z rot = — 1 they rotate clockwise. 
The parameter a determines the distribution of the stars over orbits 
with different energies whereas (3 influences the way orbits with 
different angular momenta are populated. The energy cutoff E > 
Eo will prevent stars from reaching beyond a certain maximum 
radius. The angular momentum cutoff \L z i\ > Lo on the other 
hand will impose a minimum radius that stars can reach. 



4. 3. 1 The component 's frame of reference 

These components are to be employed in the modelling of slowly 
rotating, approximately spherical galaxies using a severally sym- 
metric potential. Throughout this work, all dynamical quantities 
are expressed in spherical coordinates. These particular compo- 
nents however are axisymmetric and their characteristics will be 
most readily calculated in cylindrical coordinates. Therefore, we 
seek a transformation linking the spherical coordinates in which 
the properties of the galaxy as a whole are expressed to a cylin- 
drical coordinate system, attached to each individual component of 
this family. 

First, we introduce a cartesian reference frame (x, y, z) with 
its origin at the galaxy's centre. The components of a vector a in 
spherical coordinates (r, 9, ip) can be linked to its components in 
cartesian coordinates via the transformation 



/ O r \ / O x \ 

I a e J = A I a y I , 

\ a v ) \ a z ) 

with 

(sin 6 cos ifi sin 9 sin <p 

cos 9 cos ip cos 9 sin (p 
— sin ip cos p 




(6) 



(V) 



Suppose the direction of the rotation axis of the system we are 
describing, whether it be a small subsystem or the galaxy's main 
body, is given by the angles (#0, <po) with respect to the cartesian 
reference frame. The cartesian coordinate system (x',y' , z'), fixed 
to the component, can be obtained by first rotating the (x,y,z)- 
frame over an angle ipo around the z-axis, yielding the auxiliary 
(x", y" , z") -frame, followed by a rotation of the (x" , y" , z")- 
frame over an angle 9o around the j/"-axis : 



B 



with 

B = I cos 9o sin ipo 
\ — sin 



cos 9o cos po — sin ipo sin 60 cos <po 

•' 1 " ■ • cos ipo sin 0o sin ipo 

cos 60 



(8) 



(9) 



To facilitate the calculations as much as possible, we will 
henceforth use cylindrical coordinates (w,4>,z'). These are con- 
nected to the (x',y' , z') -frame by the transformation 



cos <f> — sin 
sin 4> cos <f> 
1 




(10) 
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One can now immediately make the connection between the spher- 
ical coordinates (r, 8, ip) in the galaxy's reference frame and the 
cylindrical coordinates (zu, (j>, z') in the component's frame of ref- 
erence using the matrix D = ABC : 



Or 
"0 




(12) 



oo •/ — oo 



4.3.2 The velocity moments 

The axisymmetric velocity moments for a counter-clockwise rotat- 
ing model {L z i ^ 0, i ro t = 1) are defined as 

OO fOC c OO 

F(E, L z i) v m v™v z , dv m dv<p d z , (13) 

with E = ip(\/zu 2 + z' 2 ) — and L z i = mv$. The moments of 
a clockwise rotating model (L z i ^ 0, i ro t = — 1) with the same 
parameters a and j3 are readily obtained since 

Ai7m,nW>,^) = Ct x pf mn (ip,vu). (14) 

The distribution function is an even function of u ro and v z i, there- 
fore only the moments with even I and n values will be non-zero. 
One can define anisotropic moments 

c + OO 



2tt 



dv 



A I 



dvtF{E,L z )vTvtt + \{\5) 



that are connected to the true moments by the relation 
1 11 

/4i, m ,2n(V>,^) = ~ B ( l + 2> n + 2 )^(i+n),m(V', W ) • (I 6 ) 

Performing some calculations the expression for the 
anisotropic velocity moments becomes: 

tin, m (1>,a>) = 2 a+2 " +(m+5)/2 ^r(a + l)r(n + l) (17) 
x^(^-S ) (Q+n+1)/2 

fm\ r(p + i + l) I L 



Z T(a + /3 + n + i + 3) V ^2^ 

i— 

/ T \ a + /3+n + i+2 

"(^-tc) 

x 2-Fi (—a — 7i — l,a + n + 2;a + /3 + n + i + 3; 

1 y/2vj*{il> - Ep) - Lq \ 

2 y/2uj 2 (tp -E ) J ' 

The true velocity moments then immediately follow : 



(18) 



•m j-)a+i + Ti+(m + 5) /2 



xr( Q + i)r(i±i)r(^ 



x(v /^3^ )t , + ((+n)/2 + l 

r(/3 + i + i) 



E/ m 



T(a + /3+ +i + 3) VV2ro 



x vV - 



_Lo_ 

a/2cs: 



i+/3+(i+ri)/2+i+2 






Figure 7. A plot of the spatial density in the meridional plane of models 
with r m ; n = 3 kpc, r max = 10 kpc, a = 4 and /3 = 4 (left), 8 (middle) 
and 12 (right). The boundary of this region is given by equation l21l hence 
all models fill the same volume. The higher the value of /3, the more the 
high-angular momentum orbits are populated. This causes the position of 
the peak density (marked by the dashed line) to slowly shift outwards. 



x 2 Ft -a ha-\ h 2; 



a + f3 + 



I + n 



v /2uj 2 (i>-E () )-L 
2^2zu 2 (^-E ) 



(19) 



4.3.3 Some commonly used moments 

4.3.3.1 Spatial density The spatial density is given by 

p{w, z') = p,o,o,ofa, z) = 



I 1 7TZU 



r(a + i)r(/3 + i) 

T(a + /3 + 3) 



XnFi(-a-l,a + 2;a + P + 3; 



Lo 
V2zc 



a + /3 + 2 



y/2vp 2 {il>-E )-Lo 
2y/2zu 2 {^-E ) 



(20) 



The density is only non-zero inside the volume defined by 
ip{r) — Eo > Lf)/(2zu 2 ). Evidently, the boundary of this volume, 
where the density becomes exactly zero, is determined by 



j 2 



2zu 



(21) 



For non-zero Lo and Eo, this volume has approximately the shape 
of a torus. If Lo = 0, the density is non-zero inside the spherical 
region defined by r ^ tp~ (Eq). All models that share the same 
values for Eo and Lo fill the same volume, independent of a or (3 
but the density distribution inside this volume does depend on a 
and /3, see figureQ 

4.3.3.2 Mean velocity For the mean rotation velocity, one finds 

Mo,i,o(^, z') 



v<p(zu,z') — 



p(tu,z') 



zu Q + p + 3\ 

2-Fi (-a - 1, a + 2; a + /3 + 4; 



V2zu 



2 F 1 ( -a - 1, a + 2; a + + 3; 

1 2^2^2 (V; _ Bo) 



(22) 



At the boundary of the component, the mean rotation velocity takes 
the value = Lo/zu = y/2(tp(r) — Eg). The parameter (3 has 
little impact on the rotation velocity profile, as shown in figure [8] 
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Figure 8. The influence of the angular momentum cut-off on the kinematics 
of models with a = 6, /3 = 2, r max = 6 kpc is demonstrated. The values 
of the inner radius are r m ; n = kpc (solid line), 0.5 kpc (dotted line), 
1.0 kpc (dashed line) and 1.5 kpc (dashed-dotted line). In the left column, a 
number of intrinsic kinematic quantities are plotted. The radius R is mea- 
sured in the equatorial plane, z 1 = 0. From top to bottom : the logarithm 
of the spatial luminosity density p(R), the radial (black) and the tangen- 
tial velocity dispersion (gray), the rotation velocity v c j > (R) and Binney's 
anisotropy parameter /3(R). In the right column, the projected kinematics 
are presented. The radius r p lies along the long axis of the models. From 
top to bottom : the logarithm of the projected luminosity density p p (r p ), 
the projected velocity dispersion cr p (r p ), the projected streaming velocity 
v p (r p ) and the Gauss-Hermite coefficients and /14 (gray). 



and figure|9| The extreme values of the rotation velocity are deter- 
mined solely by Lo and Eo, making the rotation velocity a very 
sensitive function of these parameters. 

4.3.3.3 Velocity dispersion Since I and n are completely equiv- 
alent in equation 1 1 9i . it is obvious that the velocity dispersions 
a%(vj,z') and o\i (zu, z') are always equal. This is the case for 
any model with a distribution function that depends only on E and 
L,i . One finds that 



2 ( >\ 2 1 



a + P + 3 



L 



Figure 9. The influence of the parameter (3 on the kinematics of models 
with = 6, r m ; n = 0.5 kpc and r max = 6 kpc is demonstrated. The 
values of (3 are 2 (full line), 4 (dotted line), 6 (dashed line) and 8 (dashed- 
dotted line). In the left column, a number of intrinsic kinematic quantities 
are plotted. The radius R is measured in the equatorial plane, z' = 0. From 
top to bottom : the logarithm of the spatial density p(R), the radial (black) 
and the tangential velocity dispersion (gray), the rotation velocity (R) 
and Binney's anisotropy parameter (3{R). In the right column, the projected 
kinematics are presented. The radius r p lies along the long axis of the mod- 
els. From top to bottom : the logarithm of the projected density p p (r p ), 
the projected velocity dispersion a p (r p ), the projected streaming velocity 
v p (r p ) and the Gauss-Hermite coefficients /13 (black) and hi (gray). 



2F1 -a - 2, a + 3; a + P + 4: 



2^1 — a — 1, a + 2 a + p + 3 - — 



(23) 



(a + /3 + 3)i \^/2uj 



x (yw) - Eo - 



Lq 
V2tz 



Fi ( -a - 1, a + 2; a + (3 + i + 3: 



2^/2^2 (V--B ) 



2F1 I -a- 1, a + 2; a + , 
2 / /\ 



o + 3. V 2 i^(V»-go)-£ \ 
' 2V 2 ^ 2 (*-E ) / 



(24) 
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Figures [S] and [5] also show that the anisotropy parameter 



with 

f3 p hi(zu, z 



cr ro (ro, z') 



(25) 



is almost independent of the components parameter /3 but very sen- 
sitive to the value of the angular momentum cutoff Lq. 



4.3.4 The spatial line-of-sight velocity distribution 

For the calculation of the spatial LOSVD, the DF will be directly 
integrated over the two velocity components in the plane of the 
sky. The integration over the line of sight will have to be performed 
numerically since it depends on the shape of the potential, which 
is generally only known numerically. We decompose the velocity 
vector in a component parallel to the line of sight, v p , and two com- 
ponents in the plane of the sky, Vi and V2. The polar unit vector 
can be written with respect to the vectors (e z , ei, e*2) with e 9 the 
unit vector along the line of sight and ei and e*2 unit vectors in the 
plane of the sky : 



= b\e 3 + fe 2 ei + 6362, 



(26) 



with b\ + b\ + b 2 = 1. We now rewrite the distribution function as 
a function of v p , vi and V2 '. 

F{vj, z, v^,v^,, t> z ) dv-nj dv^ dv z = F{vj, z, v p ,v\,V2) dvp dvi dv2 

= |r72 (l> - E ° - \i V l + v l + W 2)) 



irot(6i«p + &2«i + &3U2) — —) dv p dvi dv2- 



(27) 



Here and in the following, the power on the angular momentum will 
be assumed to be an integer, denoted by n. LOS VDs of components 
with non-integer values for this parameter will have to be calculated 
numerically. The velocity components in the plane of the sky, vi 
and 1)2, can be expressed in polar coordinates (v, i?) 



Vi 
V2 



V cos 1? 
v sin 



v > 0, G [0, 2tt] 



(28) 



Taking this transformation into account and introducing the nota- 
tion 

b 2 = 6cosx b ^ 0, x G [0,2tt] 

63 = frsinx, (29) 
following expression can be obtained 

F(uj,z,v,'&,Vp)dvpvdvd& = —j^ I tp - E — -{Vp + v 2 )\ 



x ( irotbiVp + bv cos J dvpvdvd-d, 

in which \ nas been dropped, since the sign and the zero-point 
of i9 are immaterial. The cosines 6 and foi obey the relation b 2 + 
b\ — 1, with b always positive. The distribution function has to be 
integrated over that part of the rectangle 



(30) 



< v < a/2(V> - E ) - vl 
s£ 1? sC 2tt 



where 



L 



iiotbiVp + bv cos 1? > 0. 



(31) 



(32) 
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Figure 10. The grey shaded area in the (v, i9)-plane marks the area over 
which the distribution function has to be integrated to calculate the spatial 
LOSVD if V > 0. The hatched rectangle is defined by < v < A and 
SC 1? $5 7r. The integration area is bounded at the right by the curve 



To relieve the notation throughout the remaining calculations, the 
following parameters will be used 

A = y/2{^-Ek)-vl, 

P = «rotOl Vp , 

ZD 

V 

cosi9 = --JT- (33) 
Ab 

The calculations will proceed along different paths depending 
on whether V ^ or V ^ 0. The area in the rectangle <3H over 
which the distribution function has to be integrated is further lim- 
ited by the positivity constraint on the angular momentum, 1321 . 
This translates in the extra constraints 



V 



bcost? 
V 



cosi? ^ 0, 
cost? < 0. 



boost? 

The spatial LOSVD can then be calculated as 

jVp 



(34) 



(A 2 



2\ot 1 2 

v ) av 



2& + ri/2 

x / ('P + 6ucosi9) n di9, 
Jo 

with the integration interval defined by <31> and 1341 . 



(35) 



4.3.4.1 V ^ In this case, the distribution function has to be 
integrated over the grey shaded area in figure Hoi Only if V < Ab 
will the second of the curves J34I intersect the rectangle <31> . 

V ^ Ab ^ The ^'-component of the angular momentum is 
positive everywhere inside the rectangle <3H . The integration over 
■& can be performed over the entire interval [0,7r], and so 



(V + bv cosi9)"di? 



(36) 
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, ■ i n 1 — n _ ( bv x " 
vViF^--,— ,1,' 



(37) 



where [n/2] is the largest integer, less or equal to n/2. Now, the 
integration over v can be tackled. This integral reduces to a Beta- 
function and the spatial LOSVD can be written in the following 
elegant form : 



2tt vj" 1 
2 n / 2 a + 

2 \ " + 1 



l( 



trotblVp 



L 



n 1-n _ b 2 (2(4>-E )-v%) 

— i \ a + Z: = 

2 2 {llotblVp _ 



(38) 



which is finite and well-behaved for all integer n and real a. 



V ^ Ab In this case, cos $o ^ and t?o £ [Si 71 "]- The sec- 
ond of the curves defined by 1341 cuts out a piece of the rectangle 
( 13 1 i and hence the distribution function has to be integrated over 
the shaded area in figurefTol: 



[ZU,Z ,v v ) = 



2<y + n/2 y j 

i=0 



(costf) l <M / (w 2 ) l/2 (.4 2 -« 2 ) a du 2 



+ 



(costf) l dtf 



Ucos f) J 



(w 2 ) l/2 (^ 2 -« 2 ) Q d^ 2 



(39) 



/>2 



/ 



The first auxiliary function fii is simply a Beta-function 

hx = A 2{a+1)+l B(a + l, l - + l). 

The second function hi is an incomplete Beta-function and can be 
rewritten in terms of a hypergeometric function : 



(40) 



ht = 



2A Z 



(JP) 

Ucostf/ 



2 + i 
x 2 Fi I -a, 1 



i+2 



L-2 i -■ f p 

2' 2'V.4&costf 



(41) 



Thus, we obtain the following expression for the spatial LOSVD 

4>(gj,z' ,v p ) = 
r(a + l) nj 2 (a+I u»vMiN) f^V 

2Q + n/2 ^ ^ ' Z^U ) r(Q + 2+ l) VP 

i=0 v 7 v 2/ 

\ ft'j. i 2ro " .2 a P" +2 1 (n 

(cost?) dtf + - , 7- ^4 p5 > - . 

^ fc!(2 + i/2) fc U&J /,„ 

We first turn our attention to the second term of 1421 . Here, only 
negative odd powers of cos # appear in the integral. The integra- 
tion can be done analytically with the aid of Gradshteyn & Ryzhik 



(1965) and yields the result : 



(cos #) 



-2(fe+l) 



(It) 



sin #0 



2k + i^a-k) j 



2(j-fc)-l 



(43) 



j=0 



As to the first term of I42t . the even and odd powers of cos 1? have 
to treated separately 



V- fn\ r(l + f) (AbV 



r(a + 2+ i) V P 



(cos tf^dtf 







f(^) 2Fl i2'^' a + 2 'VP / 



+ sin i?o 



r(i+0 / ,A6 
^ r(j) r( a + i + 2) V p 

E^^costfor 2 - 1 

i = 



(i-i) 



(44) 



Thus, we are lead to the following expression for the spatial 
LOSVD: 



f>(zu,z',v p ) 
2i9 w 



2^/2 a _|_ 



2 \ i+l 



n 1-n 6 2 (2(^-i5o)-^ 2 ) 
x 2 ^i [ — — , — - — ;a + 2: 



2' 2 



+ 2r(a + l)sint?o^ U-Bo-^J ( 



2 \ "+ 1 



irotblVp 



T \ n 



Ab cos #0 

'y^' ((l-j)/2)» 



(cos ^0)" 



!=0 



^ i + 2 \ i 



(- a ) fc (l+|) fc (i), / irot&l^-^ 



^ (|)fc(2+ f ) fe fe! ^fecos^o^^-^o)-^ / 
X2F1 (-fc,l;i-fc;(co S1 9o) 2 ) . (45) 

In the limit V —> Ab, one finds that #0 — > t and the complicated 
expression 1451 reduces to the far easier formula I38t . as it should. 
The rather daunting formula 1451 will only be well-behaved for in- 
teger values of a since then the last summation (with summation 
variable k) terminates after a finite number of terms. 

4.3.4.2 V < In this case, cos$o — — "S; > and hence 
#0 £ [0, 5 ]. The distribution function now has to be integrated over 
the shaded area in the (v, i?)-plane presented in figure fTTI Only if 
\V\ < Ab will the first of the curves 1341 intersect the rectangle 
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Figure 11. The grey shaded area in the (v, i9)-plane marks the area over 
which the distribution function has to be integrated to calculate the spatial 
LOSVD in the case V < 0. The hatched rectangle is defined by $C v ^ A 
and 1? ^ 7r. The integration area is bounded by the curve v r ~~a ■ 



\V\ ^ Ab The integration area defined by the equations <34t 
and 13 1 i vanishes, hence the spatial LOSVD is identically zero for 
all \V\ 5? Ab. 

\V\ < Ab The distribution function must be integrated over 
the grey shaded area in figure fTTI : 

d>(zo,z' ,v v ) = 



i=0 



( _!ZJ_)2 



(v 2 y /2 (A 2 -v 2 ) a dv 2 . 



(46) 



We start with the integration over the velocity in the plane of the 
sky, v, which can be performed analogously to J40> and J4H : 



{v 2 y /2 (A 2 -v 2 ) a dv 2 



A 



2(a + l)+i 



B(a + 1,- + 1) 



i + 2 \ 6 cos i9 



Thus, we obtain an expression that is almost exactly like 1421 



r(a+ ; 1) ro "^ 2(Q+1) P" 



(47) 



E 



r(i + |) /^v 



i / T(a + 2 + 



2^7" ,o„V 



(-1)' 



(/(? 



2 a +"/ 2 6 2 ^ i + 2 \ i 

i=0 



^ fci(2+i/2) fc y 1 ' 

The integrations over i? are handled analogously to <43> and J44I 
and we finally obtain 



{m,z ,v p ) = 



2i?o 

2 n/2 



— irotOlUp V - ^0 — 

a + 1 V vj J \ 1 



n l-n b 2 (2(<p-E )~v 2 p ) 

X 2-Tl , \0L-\-Z\ 7 

n / 2 \ "+ 1 

+ 2r( a + l)sintf |L- U-Eo-f 



2 \ Q+l 



T \ n 



X IrotblVp 

V vj 



_4focos #0 



4r( 0' + !) ! Vl^otbiUp 



y ^-^'(cos^o)- 2 '- 1 



i=0 v 2 



^— ' i + 2 V i 

i=0 



(-«)*(!+ 4)fc(|)fc 



iiotbiVp — ^ 



^ C|)fe(2 + |)fcfc! ^6cost? \/2(^-So)-^ 



^-fc, 1; - - fc; (costf ) 2 ) 



(49) 



This formula will only be used for integer values of a, as discussed 
in the case of 145 \ . In the limit \V\ — > Ab, $o vanishes. The spatial 
LOSVD becomes identically zero and remains zero for all \"P\ ^ 
Ab. 



5 DYNAMICAL MODEL 

The two most prominent absorption lines of the Ca II triplet were 
used as input for the modelling. The spectra were rebinned to reach 
sufficiently high S/N. The S/N was about 50 at radii less than 10 " 
and from these spectra 1 82 points were included in the fit. For larger 
radii, the S/N was about 30 and only 91 points per spectrum were 
included in the fit. 



5.1 Potential 

The observed R-band surface brightness profile was deprojected 
into a spherical luminosity density. Assuming a mass-to-light ratio, 
this lead to a mass density that was used to calculate a potential 
by means of the Poisson equation. A model with constant M/L did 
not result in a satisfactory fit. Hence, a number of models with a 
dark matter halo were considered. For all these models, the dark 
halo component was added in the form of a spherical distribution 
of mass, where the amount of matter is increasing as a power of 
the radius. Different models have different fractions of total matter 
mass to light emitting mass at a chosen radius (here 2 kpc or 0.9r e , 
this radius is only for scaling and has no physical meaning) and 
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Figure 12. In the lower panel: Mass profile for the best fitting model based 
on the Ca II triplet feature (solid line) and for the best fitting model with- 
out dark matter (dashed line). In the upper panel: the fraction of luminous 
matter relative to all matter for the best fitting models with dark matter. 

different total masses. In this way, a grid of models was created, 
and the smallest value for % 2 was taken to be indicative for the best 
matching potential. 

The best results were obtained with a potential where the spa- 
tial density of the total amount of matter was 2.5 times the spatial 
density of the amount of luminous matter within 2 kpc. This scal- 
ing implies that the total mass within lr e is 1.6 x 1O 11 M0, 70 % 
of which is luminous matter. The total mass at 2r e is 5 x 1O 11 M0, 
and 37 % of this is luminous matter (see also upper panel in figure 
1 1 2t . Fi sure ITzl (lower panel) shows the total mass profile in solid 
line in comparison with the luminous matter mass profile in dashed 
line. 

It was not necessary to include a central black hole in the 
model to obtain a good fit for the central spectra. 

5.2 Dynamical model 

The model that best fits the spectra around the two most promi- 
nent Ca II triplet lines, at 8542 A and 8662 A, is shown in figure 
om these plots, no obvious shortcoming of the model can be 
indicated. 

Figure fT4l shows some derived quantities of the model. The 
upper left panel shows the projected density on major and minor 
axis. As can be seen, the flattening of the galaxy is well repro- 
duced. The right panels show calculated projected mean velocity 
and velocity dispersion profiles (solid lines), together with uncer- 
tainties on these quantities as derived from the model. The values 
and error estimates for the projected mean velocity (v p ) and veloc- 
ity dispersion a p , presented in 12.41 are plotted on top in dots. In 
figure rT4l also the anisotropy parameters (3$ = 1 — (o^/cr^) (in 
solid line) and (3$ = 1 — (el/of) (in dotted line) are displayed in 
the lower left panel. The parameter /3^, Binney's anisotropy param- 
eter, is close to zero for the central parts, indicating an isotropic DF 
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Figure 14. Some derived moments from the model. Left column: projected 
density on major and minor axis and the anisotropy parameters Pa and (3g. 
Right column: projected mean velocity and projected velocity dispersion. 
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Figure 15. Template mixes for the best model based on the two strongest 
lines of the Ca II triplet. 



there, becomes slightly posivite for intermediate radii and slightly 
negative for radii larger than 2 kpc, indicating that the model is 
tangentially anisotropic there. 

A mix of G5 V and K4III stars gives the best fit to the spectra in 
this wavelength range. The value for % 2 obtained for the fit is 2267, 
the second best model has a \ 2 of 2383. The resulting template mix 
of this best model can be found figure lT5l The relative densities are 
calculated from the spatial number densities of the model. 

More models, with other template mixes, than shown here 
were calculated. The results mainly confirm the trends seen in fig- 
ure^] In the model with K4III and G2V stars e.g., the contribution 
of G2V stars to the density is larger than the G5 V stars in the K4III- 
G5V panel, but the same radial behaviour occurs. 



6 DISCUSSION 

If a model is strongly tangential in the outer regions, most of the 
motion occurs along the line of sight and this may increase the pro- 
jected velocity dispersion. Hence, the influence of dark matter may 
be underestimated. This is the strongest argument for including the 
fourth order shape parameter of the observed LOSVD as input pa- 
rameter in a dynamical model, because it is believed to indicate 
whether a lot of orbits are observed near their tangent point or not. 
For the calculation of this dynamical model, we did not rely on the 
extracted LOSVD nor the derived kinematic parameters, but on the 
observed spectra. 

The dynamical model is slightly tangentially anisotropic from 
roughly lr e on, as can be seen in figure fT4l 

Tangential anisotropy in the outer regi ons of the model is als o 
seen in other models for e.g. NGC 3377 jGebhardt etai]l200fj|) . 
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Figure 13. Data and model around the two strongest Ca II triplet lines. The regions included in the fit are located between the intervals indicated by dashed 
lines and with centre around 8620 A and 8750 A. The projected radius of the positions where the spectra are taken are indicated in the panels. 
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Figure 16. LOSVDs as reconstructed from the kinematic parameters pre- 
sented in l2.4l (dashed line) and calculated from the model (solid line). Also 
eiTor bars are calculated from the model. The radii at which the LOSVDs 
are taken are indicated in arcsec in the panels. 



Figure 17. LOSVDs calculated from the model: for the composed DF in full 
lines, for the G5V population in dotted lines and for the K4III population 
in dashed lines. The different radii at which the LOSVDs are taken are 
indicated in arcsec in the panels. 



NGC 1023 bower et afcOOlh and also in the sample modeled by 
Kron awitter et alJfcOOOl) . 

Estimating kinematic parameters through a dynamical model 
clearly reduces the amount of scatter in the kinematic profile, as 
is illustrated in figure fl4l fright panels). This is easily understood 
since the model induces additional dynamical constraints on the 
parameters and as a result the values of the parameters taken at dif- 
ferent radii are not independent. When the parameters are inferred 
directly from the observed spectra, there is no connection between 
the parameters at different radii. 

Figure fTrjl presents LOSVDs as derived from the observations 
and the model. The profiles calculated from the kinematic param- 
eters presented in 12.41 are shown in dashed lines, the solid lines 
present the LOSVDs with error bars calculated from the model. It 
is clear that the LOSVDs calculated from the kinematic parameters 
all have an approximately Gaussian shape. This is of course due 
to the parametrization that is used to present them. The LOSVDs 
calculated from the model show a large variety of shapes. Some 
of these profiles cannot be reproduced by a Gauss-Hermite series 
truncated at fourth order (or even sixth order) and hence a perfect 
agreement with the LOSVDs calculated from the kinematic param- 
eters cannot be expected. 

The shapes of the LOSVDs calculated from the model can 
be better understood when looking at figure fT71 In this figure the 
LOSVDs for the different stellar populations, together with the 
composed LOSVD (solid lines) are shown. The relative height of 
the LOSVDs for the two populations is in agreement with the in- 
formation in figure [131 A full discussion of the role of the stellar 
populations will be deferred to a subsequent paper, where a fit will 
be presented that takes into account 2 wavelength regions. 

Instead of comparing full profiles, one could argue in favor 
of comparing their most relevant parameters. Unfortunately, this 
is not a trivial exercise. The mathematical elegane of the Hermite 
expansion lies in the fact that it forms an orthogonal series, with 



coefficients that can be calculated by integrals that follow from the 
theory. However, this is almost never the way the Hermite coef- 
ficients are calculated in practice: they are generally determined 
with some kind of minimization after performing a convolution. 
This causes the loss of the orthogonality property. This becomes 
clear in an exercise where a stellar template spectrum convolved 
with a LOSVD with parameters — 20km/s < (v) < 20km/s, 
a = 300km/s , h$ — 0, hi — 0, he — and he = 0.1 is analysed 
in the same way as real galactic spectra and kinematic parameters 
up to hi are determined. Since the template spectrum used to cre- 
ate this synthetic spectrum is equal to the one used to determine the 
parameters again, there is no influence of continuum subtraction or 
template mismatch on the results. The synthetic galaxy spectrum 
has a S/N of 130 in the centre and the light profile follows a de 
Vaucouleurs that resembles the light profile of NGC 3258. The re- 
sult of this simulation is shown in figure fT8l The values retrieved 
for (v) and /13 correspond to the original values. On the other hand, 
the values for a and hi are clearly higher than the values of the 
input LOSVDs. 

This indicates that the kinematic parameters as derived from a 
Gauss-Hermite parametrization are in practice not independent, in 
the sense that power in the higher orders than considered in the ex- 
pansion may contaminate lower order parameters. In the case of hi 
this can give rise to wrong conclusions concerning the anisotropy 
in the stellar motion. Naively, one could conclude that it is better 
to go to higher order than hi in the parametrization. Unfortunately, 
the signal-to-noise ratio in outer regions of a galaxy is hardly suf- 
ficient to derive reliable values for hi, so it does not make sense to 
try to determine higher order parameters. The useful information in 
LOSVDs derived from observations is thus very limited, due to the 
rather low signal-to-noise of the observations. Therefore, though it 
is a good idea to use a parametrization with a limited number of 
parameters to represent these LOSVD, there is, as a consequence, 
also only a limited number of shapes for LOSVDs possible, as ex- 
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Figure 18. Result from a simulation where (v), a, hg and /14 are deter- 
mined from spectra that are a convolution of LOSVDs with — 20km/s < 
(v) < 20km/s, cr = 300km/s, h 3 = h 4 = h 5 = and h 6 = 0.1. 

traded from these observations. This on top of the fact that these 
parameters are not independent, which means that caution should 
be taken when interpreting them. On the other hand, LOSVDs cal- 
culated from dynamical models can have a much larger variety of 
shapes. It would be unfair to reject model LOSVDs because they 
do not match the ones that are derived from observations, simply 
because the information in LOSVDs derived from spectra is very 
limited, while the dynamical models offer much more freedom. 



7 CONCLUSIONS 

In this paper we use a method to analyse spectra from elliptical 
galaxies in order to retrieve dynamical information and to perform a 
stellar population synthesis at the same time. The method originates 
from the field of dynamical modelling and is presented by DD98. 
The main idea is that the different stellar populations that contribute 
to the integrated galaxy spectrum do not necessarily share the same 
kinematic characteristics. Hence, in terms of dynamical modelling, 
they should be given different distribution functions. This paper 
reports for the first time on the application of our modelling tech- 
nique to data. The attention is focused on the presentation of the 
method and on establishing a mass estimate for the galaxy under 
study. A second paper will discuss the possibilities for stellar popu- 
lation studies and will go deeper into the distribution functions for 
separate stellar templates. 

Modelling the spectra directly yields an estimate for the to- 
tal amount of mass, and an analysis of the internal dynamics of 
the galaxy. This technique, using galaxy spectra as input data, also 
differs from more traditional modelling strategies that rely on kine- 
matic parameters that are derived from the spectra in the sense that 
a two step process is turned into a one step process, thereby coming 
closer to the astrophysical issues. 

Photometric and spectroscopic observations for NGC 3258 are 
presented. This galaxy has almost perfect elliptical isophotes, and 
has an effective radius of 12.95" (2.26 kpc). Images in B and R 
show a slightly off-centre nucleus. There seems to be a slight red- 



dening of the nucleus. The photometry also reveals the presence of 
a central dust disk (0.17 kpc - 0.26 kpc). The R band photometry 
was used to derive a spherical potential for the luminous matter. 

The photometry, together with galaxy spectra containing dy- 
namical information up to 1.3r e are combined into a dynami- 
cal model. The best dynamical model for NGC 3258 is the one 
with a dark matter halo that contributes 30% of the total mass 
at lr e and 63% of the total mass at 2r e . From 1.45r e outwards, 
the model is dominated by dark matter. The total mass within 
lr e is 1.6 x 10 11 A'/q and (M/L) R w 65, (M/L) R w 30, 
whereas the total mass extrapolating to 2r e is 5 x 10 11 Mq and 
(M/L) R w 95,(M /L) B « 43.5. In the region without dark mat- 
ter (M/L)r w 14 ((M/L) B ~ 6.4). The modelling did not re- 
quire to include a central black hole. 

The anisotropy parameter (see figure 1 14i indicates that the 
model is isotropic in the centre, is radially anisotropic between 
0.2 and 2 kpc (i.e. within lr e ) and becomes slightly tangentially 
anisotropic further out. 

A simulation has shown that for the analysis of LOSVDs de- 
rived from observations by means of a truncated Gauss-Hermite se- 
ries, power in higher orders than considered in the expansion may 
contaminate lower order parameters. This is a wortysome result, 
since hi profiles derived from observations tend to be interpreted in 
terms of anisotropy of the galaxy. If this parameter is contaminated 
by power in higher orders by some form of aliasing induced by 
the fitting process, this may lead to an erroneous interpretation or 
at least may make comparison of /14's obtained by different meth- 
ods problematic. Moreover, due to the freedom in the dynamical 
model, one cannot expect that LOSVDs calculated from a dynami- 
cal model as presented in this paper can be approximated by Gauss- 
Hermite series that are truncated at a low order. 
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